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Energetics and conductance in jellium modelled nanowires 
are investigated using the local-density-functional-based shell 
correction method. In analogy with studies of other finite-size 
fermion systems , e.g., simple-metal clusters or ^He clusters, 
we find that the energetics of the wire as a function of its 
radius (transverse reduced dimension) leads to formation of 
self-selecting magic wire configurations (MWC's, i.e., discrete 
sequence of wire radii with enhanced stability), originating 
from quantization of the electronic spectrum, namely forma- 
tion of subbands which are the analogs of electronic shells in 
clusters. These variations in the energy result in oscillations 
in the force required to affect a transition from one MWC of 
the nanowire to another, and are correlated directly with step- 
wise variations of the quantized conductance of the nanowire 
in units of 2e^ /h. 



I. INTRODUCTION 

Identification and understanding of tlie physical origins 
and systematics underlying the variations of materials 
properties with size, form of aggregation, and dimension- 
ality are some of the main challenges in modern materi- 
als research, of ever increasing importance in the face of 
the accelerated trend toward miniaturization of electronic 
and mechanical devices. While for over two decades stud- 
ies of size-evolutionary patterns oimaterials have focused 
on atomic and molecular clustersou in beams or embed- 
ded in inert matrices, more recent efforts concentrated 
on preparation, characterization, and understanding of 
finite solid-sate structures. These include jjiapometer- 
scalc metal and semiconductor nanocryatals,QQ surface- 
supported structures and quantum dots|3 and nanoscale 
junctions or wires. ^~^^ 

Interestingly, it has emerged that concepts and 
methodologies developed in the context of isolated gas- 
phase clusters and atomic nuclei are often most useful 
for investigations of finite-size solid-state structures. In 
particular, it has been shown most recentlyOEZl through 
first-principles molecular dynamics simulations that as 
metallic (sodium) nanowires are stretched to just a few 
atoms in diameter, the reduced dimensions, increased 
surface-to-volume ratio, and impoverished atomic envi- 
ronment, lead to formation of structures, made of the 
metal atoms in the neck, which can be described in 
terms of those observed in small gas-phase sodium clus- 
ters; hence they were termed as supported cluster- derived 
structures (cds). The above prediction of the occurrence 
of "magic-number" cds's in nanowires, due to charac- 
teristics of electronic cohesion and atomic bonding in 



such structures of reduced dimensions, are directly corre- 
lated with the energetics of metal clusters, where magic- 
number sequences of clusters sizes, shapes and structural 
motifs due to electronic and/or geome^tcie-shell effects, 
have been long predicted and observed.BBua 

Furthermore, these results lead one directly to con- 
clude that other properties of nanowires, derived from 
their energetics, may be described using methodologies 
developed previously in the context of clusters. In- 
deed, in this paper, we show that certain aspects of 
the mechanical response (i.e., elongation force) and elec- 
tronic transport (e.g., quantized conductance) in metal- 
lic nanowires can be analyzed using the local-density- 
approximation (LDA) -based shell correction method 
(SCM), develqaed and applied previously in studies of 
metal clusters.oO Specifically, we show that in a jellium- 
modelled, volume-conserving nanowire, variations of the 
total energy (particularly terms associated with elec- 
tronic subband corrections) upon elongation of the wire 
lead to self- selection of a sequence of stable "magic" wire 
configurations (MWC's, specified in our model by a se- 
quence of the wire's radii), with the force required to 
elongate the wire from one configuration to the next ex- 
hibiting an oscillatory behavior. Moreover, we show that 
due to the quantized nature of electronic states in such 
wires, the electronic conductance varies in a quantized 
step-wise manner (in units of the conductance quantum 
5o — 2e^//i), correlated with the transitions between 
MWC's and the above-mentioned force oscillations. 

Prior to introducing the model studied in this paper, it 
is appropriate to briefly review certain previous theoreti- 
cal and experimental investigations, which form the back- 
ground and motivation for this study. Atomistic descrip- 
tions, based on realistic interatomic interactions, and/or 
first-principles modelling and simulations played an es- 
sential role in discovering the formation of nanowires, 
and in predicting and elucidating the microscopic mech- 
anisms underlying their mechanical, spectral, electronic 
and transport properties. 

Formation and mechanical properties of interfa- 
cial junctions (in the form of crystalline nanowires) 
have been predicted through early molecular-dynamics 
simulations^ where the materials (gold) were modelled 
using semiempirical embedded-atom potentials. In these 
studies it has been shown that separation of the con- 
tact between materials leads to generation of a connective 
junction which elongates and narrows through a sequence 
of structural instabilities; at the early stages, elongation 
of the junction involves multiple slip events, while at the 
later stages, when the lateral dimension of the wire necks 



down to a diameter of about 15 A, further elongation in- 
volves a succession of stress accumulation and fast relief 
stages associated with a sequence of order-disorder fttaVn 
tural transformations localized to the neck region .HO'ES 
These structural evolution patterns have been shown 
through the simulations to be portrayed in oscillations 
of the force required to elongate the wire, with a period 
approximately equal to the interlayer spacing. In addi- 
tion, the "sawtoothed" character of the predicted force 
oscillations [see Fig. 3(b) in Ref. || and Fig. 3 in Ref. |2|] 
reflects the stress accumulation and relief stages of the 
elongation mechanism. Moreover, the critical resplyed 
yield stress of gold nanowires has been predictedaO to 
be ~ 4GPa, which is over an order of magnitude larger 
than that of the bulk, and is comparable to the theoreti- 
cal value for Au (1.5 GPa) in the absence of dislocations. 

These predictions, as,rwell as anticipated electronic 
conductance propertiespfl have been corroborated in a 
number of experupents . UpBgi scanning tunnclipg and 
force microscopy JsliHdiMSEjGfl break j unctions Jlj and 
pin-plate techniquesEjo at ambient environments, as 
well as under ultrahigh vacuum and/or cryogenic con- 
ditions. Particularly, pertinent to our current study are 
experimental observations of the oscillatory behavior of 
the elongation forces and the correlations between the 
changes in the conductance and the force oscillations; see 
especially the simultaneous measurements of force and 
conductance in gold nanowires in Ref. |l^, where in addi- 
tion the predicted "ideal" value of the critical yield stress 
has also been measured (see also Ref. 24 ) . 

The jellium-based model introduced in this paper, 
which by construction is devoid of atomic crystallo- 
graphic structure, does not address issues pertaining to 
nanowire formation methods, atomistic configurations, 
and mechnanical response modes [e.g., plastic deforma- 
tion mechanisms, interplanar slip, ordering and disorder- 
ing mechanisms (see detailed descriptions in Refs. |qj2^ 
and p3, and a discussion of conductance dips in Ref. 12), 
defects, mechanichal reversibility,E3L3 and-xoughening of 
the wires's morphology during elongationc3| , nor does it 
consider the effects of the above on the electron spec- 
trum, transport properties, and dynamic£3 Neverthe- 
less, as shown below, the model offers a useful framework 
for linking investigations of solid-state structures of re- 
duced dimensions (e.g., nanowires) with methodologies 
developed in cluster physics, as well as highlighting cer- 
tain nanowire phenomena of mesoscopic origins and their 
analogies to clusters. 



II. THE JELLIUM MODEL FOR METALLIC 

NANOWIRES: THEORETICAL METHOD AND 

RESULTS 

Consider a cylindrical jellium wire of length L, hav- 
ing a positive background with a circular cross section 
of radius R <S^ L. For simplicity, we restrict ourselves 



here to this symmetry of the wire cross section. Vari- 
ations in the shape of the nanowire cross section sersie 
to affect the degeneracies of the electronic spectrurnc3'E£l 
without affecting our general conclusions. We also do 
not include here variations of the wire's shape along its 
axis. Adiabatic variation of the wire's axial shape intro- 
duces a certain amount of smearing of the conductance 
steps through tunnellina,-denending on the axial radius 
of curvature of the wire.EJ'Eilc3 Both the cross-sectional 
and axial shape of the wire can be included in our model 
in a rather straightforward manner. 

The principal idea of the SCM is the separation of the 
total LDA energy EriR) aMtBM 



Et{R) = E{R) + AE,Y,{R) , 



(1) 



where E{R) varies smoothly as a function of the sys- 
tem size, and AEsh{R) is an oscillatory term arising 
from the discrete quantized nature of the electronic lev- 
els. AEsh^) is usually,|-Galled a shell correction in 
the nuclearE2l and clusterac3 literature; we continue to 
use here the same terminology with the understand- 
ing that the electronic levels in the nanowire form sub- 
bands, which are the analog of electronic shells in clus- 
ters, where furthermore the size of the system is usu- 
ally given by specifying the number of atoms iV. The 
SCM method, which has been shown tn_4deld results 
in excellent agreement with|-ratperimenta3'EJ'Ej and self- 
consistent LDA calculationaJ'CJ for a number of cluster 
systems, is equivalent to a Harris functional (-Eharris) 
approximation to the Kohn-Sham LDA with the in- 
put density obtained through variational minimization 
of an extended Thomas-Fermi (ETF) energy functional, 
EetfIpI (with the kinetic energy, Tetf[p]' given to 4th 
order gradients and the potential, Vetf, including the 
Hartree repulsion and exchange-correlation and positive- 
background attractions as in LDA). The smooth con- 
tribution in Eq. J^ is identified with E-etf[p\- The 
optimized densityEj p at a given radius R is obtained 
under the normalization condition (charge neutrality) 

2T:Jp{r)rdr = P^l^{R), where p^^\r) = 3i?V(4r3) is 
the linear positive background density. Using the opti- 
mized p, one solves for the eigenvalues e^ of the Hamilto- 
nian H — — (?i^/2m)V^ + Vetf[p], and the shell correc- 
tion is given by 

AEsh = Ehan-is [p] — -E^ETF [p\ 
occ ^ 

= E ^* - / ?(i-)^ETF [p{r)]dr " Tetf [p\ , (2) 

i—l 

where the summation extends over occupied levels, here 
the dependence of all quantities on the pertinent size vari- 
able (i.e., the radius of the wire R) is not shown explicitly. 
Additionally, the index i can be both discrete and con- 
tinuous, and in the latter case the summation is replaced 
by an integral. 

Following the above procedure with a uniform back- 
ground density of sodium (r^ ~ 4 a.u.), a typical po- 



tential Vetf('") for R = 12.7 a.u., where r is the radial 
coordinate in the transverse plane, is shown in Fig. 1, 
along with the transverse eigenvalues ?„„ and the Fermi 
level; to simplify the calculations of the electronic spec- 
trum, we have assumed (as noted above) R ^ L, which 
allows us to express the subband electronic spectrum as 
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where kz is the electron wave number along the axis of 
the wire (z). 

As indicated earlier, taking the wire to be charge neu- 
tral, the electronic linear density, p)^ (R), must equal the 

linear positive background density, pj^ (R)- The chemi- 
cal potential (at T = the Fermi energy ep) for a wire 
of radius R is determined by setting the expression for 
the electronic linear density derived from the subband 
spectra equal to p)^ (R), i-c. 
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where the factor of 2 on the left is due to the spin de- 
generacy. The summand defines the Fermi wave vec- 
tor for each subband, kp.mn- The resulting variation 
of epiR) versus R is displayed in Fig. 2(a), showing 
cusps for values of the radius where a new subband 
drops below the Fermi level as R increases (or conversely 
as a subband moves above the Fermi level as R de- 
creases upon elongation of the wire). Using the Lan- 
dauer expression for the conductance G in the limit of 
no mode mixing and assuming unit transmission coef- 
ficients, G{R) = goY.n,m'^i'^FiR) - e«m(i?)], wherc 6 
is the Heaviside step function. The conductance of the 
nanowire, shown in Fig. 2(b), exhibits quantized step- 
wise behavior, with the step-rises coinciding with the lo- 
cations of the cusps in ep{R), and the height sequence of 
the steps is Igo, 2go, 2go, Igo, ..., reflecting the |eircular 
symmetry of the cylindrical pares' cross sections,a as ob- 
served for sodium nanowires.Ej Solving for ep{R) [see Eq. 
(0)], the expression for the sum on the right-hand-side of 
Eq. (H) can be written as 



occ „ occ ^kf 



dkzC'am.y'^z'i ^) 



— YypiR) + 2?„,„(i?)]^^[ef (i?) - ?-„„,(i?)] , (5) 

which allows one to evaluate Ai?sh [Eq- (||)] for each wire 
radius R. Since the expression in Eq. (^) gives the energy 
per unit length, we also calculate E'etf, 2etf7 and the 
volume integral in the second line of Eq. (0) for cylin- 
drical volumes of unit height. To convert to energies per 
unit volume [denoted as exiR), £{R), and AeghiR)] a-^ 
energies are further divided by the wire's cross-sectional 



area, wR^. The smooth contribution and the shell correc- 
tion to the wire's energy are shown respectively in Fig. 
3(a) and Fig. 3(b). The smooth contribution decreases 
slowly towards the bulk value (—2.25 eV per atomcS). 
On the other hand, the shell corrections are much smaller 
in magnitude and exhibit an oscillatory behavior. This 
oscillatory behavior remains visible in the total energy 
[Fig. 3(c)] with the local energy minima occurring for 
values i?niin corresponding to conductance plateaus. The 
sequence of i?niin values defines the MWC's, that is a 
sequence of wire configurations of enhanced stability. 

From the expressions for the total energy of the wire 
[i.e., rier(-R), where fi = ttR^L is the volume of the wire] 
and the smooth and shell (subband) contributions to it, 
we can calculate the "elongation force" (EF), 



Ft{R) 



dineriR)] 

dZ 



-n 



r de{R) 
{ dL 



d[Ae,^{R)] \ 
dL J 



F{R) + AF,h(i?) 



(6) 



Using the volume conservation, i.e., d{TiR^L) = 
0, these forces can be written as Ft{R) — 
iTTR^/2)dsT{R)/dR, F{R) = {TrR^/2)de{R)/dR, and 
AFshiR) = {nR^/2)d[AeshiR)]/dR. F{R) and Ft{R) 
are shown in Fig. 3(d,e). The oscillations in the force re- 
sulting from the shell-correction contributions dominate. 
In all cases, the radii corresponding to zeroes of the force 
situated on the left of the force maxima coincide with the 
minima in the potential energy curve of the wire, corre- 
sponding to the MWC's. Consequently, these forces may 
be interpreted as guiding the self-evolution of the wire to- 
ward the MWC's. Also, all the local maxima in the force 
occur at the locations of step-rises in the conductance [re- 
produced in Fig. 3(f)], signifying the sequential decrease 
in the number of subbands below the Fermi level (con- 
ducting channels) as the wire narrows (i.e., as it is being 
elongated). Finally the magnitude of the total forces is 
comparable to the measured ones (i.e., in the nN range). 



III. CONCLUSIONS AND DISCUSSION 

We investigated energetics, conductance, and meso- 
scopic forces in a jellium modelled nanowire (sodium) 
using the local-density-functional-based shell correction 
method. The results shown above, particularly, the os- 
cillations in the total energy of the wire as a function 
of its radius (and consequently the oscillations in the 
EF), the corresponding discrete sequence of magic wire 
configurations, and the direct correlation between these 
oscillations and the step-wise quantized conductance of 
the nanowires, originate from quantization of the elec- 
tronic states (i.e., formation of subbands) due to the re- 
duced lateral (transverse) dimension of the nanowires. 
In fact such oscillatory behavior, as well as the appear- 
ance of "magic numbers" and "magic configurations" of 
enhanced stability, are a general characteristic of finite- 
size fermionic systems and are in direct analogy with 



those foupd in simple-metaLj;lusters (as well as in '^He 
clusterstj and atomic Buciej — ), where electronic shell ef- 
fects on thfi energeticsu'oEaO (and most recently shape 
dynamicsLj of jcUium modelled clusters driven by forces 
obtained from shell-corrected energetics) have been stud- 
ied for over a decade. 

While these calculations provide a useful and instruc- 
tive framework, we remark that they are not a substi- 
tute for theories where the atoniistic, jiaiuLe_and spe- 
cific atomic arrangements are included3^E3oou in eval- 
uation of the energetics (and dynamics) of these sys- 
tems (see in particular Refs. p5l27, where first-principles 
molecular-dynamics simulations of electronic spectra, ge- 
ometrical structure, atomic dynamics, electronic trans- 
port and fluctuations in sodium nanowires have been dis- 
cussed). 

Indeed, the atonfistic structural characteristics of 
nanowires (including the occurrence—of cluster-derived 
structures of particular geometriesESEj) , which may he, 
observed through the use of high resolution microscopy ,E3 
influence the electronic spectrum and transport charac- 
teristics, as well as the energetics of nanowires and their 
mechanical properties and response mechanisms. In par- 
ticular, the mechanical response of materials involves 
structural changes through displacement and discrete re- 
arrangenent of the atoms. The mechanisms, pathways, 
and rates of such structural transformations are depen- 
dent on the arrangements and coordinations of atoms, 
the magnitude of structural transformation barriers, and 
the local shape of the wire, as well as possible dependency 
on the history of the material and the conditions of the 
experiment (i.e., fast versus slow extensions). Further ev- 
idence for the discrete atomistic nature of the structural 
transformations is provided by the shape of the force vari- 
ations (compare the calculated Fig. 3(b) in Ref. S and 
Fig. 3 in Ref. p2| with the measurements shown in Figs. 
1 and 2 in Ref7]l9|) , and the interlayer spacing period of 
the force oscillations when the wire narrows. While such 
issues are not addressed by our model, the mesoscopic 
(in a sense universal) phenomena described by it are of 
interest, and may guide further research in the area of 
finite-size systems in the nanoscale regime. 
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FIG. 3. (a-c): The smooth (a) and shell-correction 

(b) contributions to the total energy (c) per unit vol- 
ume of the jellium-modelled sodium nanowire (in units of 
u = 10"'' eV/a.u.^), plotted versus the radius of the wire 
(in a.u.). Note the smaller magnitude of the shell corrections 
relative to the smooth contribution, (d-e): The smooth con- 
tribution (d) to the total force and the total force (e), plotted 
in units of nN versus the wire's radius. In (e), the zeroes 
of the force to the left of the force max;ima occur at radii 
corresponding to the local minima of the energy of the wire 
(c). In (f), we reproduce the conductance of the wire (in 
units of po = 2e^//i), plotted versus R. Interestingly, calcula- 
tions of the conductance for the MWC's (i.e., the wire radii 
corresponding to the lft«s*ions of the step-rises) through the 
Sharvin-Weyl formula ,li3|'E3 corrected for the finite height of 
the confining potentialli3 (see lower panel of Fig. 1), namely 
G — go{'KS/\p — aP/Xp) where 5* and P are the area and 
perimeter of the wire's cross section and \f is the Fermi wave- 
length (Af = 12.91 a.u. for Na) with a = 0.1 (see Ref. ^, 
yield results which approximate well the conductance values 
(i.e., the values at the bottom of the step-rises) shown in (f). 



FIG. 1. Lower panel: The Vetf{j) potential for a sodium 
wire with a uniform jellium background of radius R = 12.7 
a.u., plotted versus the transverse radial distance from the 
center of the wire, along with the locations of the bottoms 
of the subbands (namely the transverse eigenvalues £„„; n 
is the number of nodes in the radial direction plus one, and 
m is the azimuthal quantum number of the angular momen- 
tum). The Fermi level is denoted by a dashed line. Top panel: 
The jellium background volume density (dashed line) and the 
electronic volume density p{r) (solid line, exhibiting a char- 
acteristic spillout) normalized to bulk values are shown. 



FIG. 2. Variation of the Fermi energy ep [shown in 
(a)] and of the conductance G (shown in (b) in units of 
go = 2e^/h), plotted versus the radius R, for a sodium 
nanowire. Note the coincidence of the cusps in ep with the 
step-rises of the conductance. The heights of the steps in G 
reflect the subband degeneracies due to the circular shape of 
the wire's cross section. 
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